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First principles calculation of the phonons modes in the hexagonal YMn03 
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The lattice dynamics of the YMnOs magneto-electric compound has been investigated using den- 
sity functional calculations, both in the ferroelectric and the paraelectric phases. The coherence 
between the computed and experimental data is very good in the low temperature phase. Using 
group theory, modes continuity and our calculations we were able to show that the phonons modes 
observed by Raman scattering at 1200K are only compatible with the ferroelectric P^cm space 
group, thus supporting the idea of a ferroelectric to paraelectric phase transition at higher temper- 
ature. Finally we proposed a candidate for the phonon part of the observed electro-magnon. This 
mode, inactive both in Raman scattering and in Infra- Red, was shown to strongly couple to the 
Mn-Mn magnetic interactions. 
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I. INTRODUCTION 

Materials presenting magneto-electric coupling exhibit 
magnetic properties coupled to the electric properties, 
such as polarization or dielectric constant. These mate- 
rials have attracted a lot of attention over the last years 
since the magneto-electric coupling allows a possible con- 
trol of the magnetic properties by an electric field and 
over electric properties using a magnetic field. 

Unfortunately, the microscopic origin of the coupling 
between the magnetic and electric order parameters is 
still ill known. The knowledge of phonons spectra can 
however bring help on understanding this coupling. In- 
deed, not only the phonons modes are strongly related to 
the existence and amplitude of a spontaneous polariza- 
tion, but in addition strong spin-phonons coupling occur 
in multiferroic material. This coupling can even be strong 
enough in order to result in hybrid excitations built from 
the mixing between phonons and spin-waves^. The exis- 
tence of such hybrid modes, called electromagnons, were 
recently discovered in RMnOs orthorhombic manganites 2 
using optical measurements. More recently such excita- 
tions were also found by inelastic neutrons scattering in 
the hexagonal manganite YMnOr^. 

Hexagonal YMnOs is a layered compound where the 
manganese ions are located in triangle-based bipyramids. 
These bipyramids are arranged in planes parallel to the 
(a, b) direction so that the manganese ions form a dis- 
torted triangular lattice (see figure [1]). They share an 
oxygen atom in the (a, b) planes. The yttrium atom 
is located in between the bipyramids layers. YMnOs 
is paraelectric at high temperature with space group 
P6s/mmc. Under a critical temperature it is ferroelec- 
tric, with space group P6scm ; the polarization is aligned 
along the c direction. The temperature of this paraelec- 
tric (PE) to ferroelectric (FE) phase transition is however 
still under debate. Indeed, usually assumed to be around 
~ 920 if- 5 , it was recently proposed to occur at a higher 
temperature, c± 1050 — 1100 K^. A very recent neu- 
trons diffraction investigation 8 suggests the existence of 
an iso-symmetric phase transition at ~ 920 K, resulting 
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Figure 1: (a) representation of the YMnOs crystal structure, 
(b) view along the c direction of the MnOs bipyramids layers. 



in a strong lowering of the polarization amplitude, while 
the non-centrosymmetric to centrosymmetric transition 
occurs at 1258 ± AK only. At much lower temperature 
(Tat — 75 K) YMnOs undergoes a paramagnetic (PM) to 
antiferromagnetic (AFM) transition. The Mn 3+ ions are 
in a high spin state S = 2. Neutrons diffraction exper- 
iments witness an in-plane spins orientation associated 
with a 120° ordering, with the spins perpendicular to 
the Mn-0 bonds^^. The ferroelectric and antiferromag- 
netic order parameters are coupled and a giant atomic 
motion was revealed at the Neel temperature—. The 
magnetic transition is however an iso-structural transi- 
tion and the space group remains P6scm in the AFM 
phase 10 . The magnetic space group is more subject to 
caution. Indeed, while most authors favor an identical 
P63C171 magnetic and structural space group2*i£, the ab- 
sence of a linear magneto-electric coupling and the exis- 
tence of a ferromagnetic component recently observed^i^ 
are incompatible with the P6scm magnetic group. Group 
analysis shows that these experimental data can only be 
accounted for in the P6' 3 magnetic group^. 

Phonons spectra are often used to help understand- 
ing both the ferroelectric/paraelectric transitions and the 
magneto-electric coupling. While the YMnOs phonons 



spectrum has been experimentally studied both by Ra- 
man^— and Infra-RecU 6 measurements, there is still 
many opened questions concerning not only modes assig- 
nations but also the fact that several of the symmetry 
allowed modes are not experimentally observed. First 
principle calculations not only could help understanding 
the experimental observation but also could yield more 
light on the recently observed electromagnon 3,4 . The 
present paper thus presents first principle calculations 
of the YMn03 phonons spectrum in the different space 
groups associated with the FE and PE phases, that is 
PQ^cm and P6s/mmc. The next section will be de- 
voted to the methodological details and preliminary re- 
sults (geometry optimization, Born effective charges, po- 
larization, dielectric tensor, etc..) and section III will 
describe and discuss phonons results. 



II. METHOD AND PRELIMINARY RESULTS 

We performed geometry optimization, spontaneous po- 
larization and phonons frequencies calculations within 
the Density Functional Theory (DFT) as implemented in 
the Crystal09 package 17 . Since the manganese 3d shells 
are strongly correlated we used hybrid functionals in or- 
der to better take into account the self-interaction cancel- 
lation. Three different functionals were used in this work, 
namely the classical B3LYP functional^, the recently 
developed B1WC— functional for ferroelectric systems 
and a B1PW— functional corresponding to the B1WC 
one with the Perdew-Wang exchange functional 20 . Only 
the B1PW results are presented in the paper, the re- 
sults obtained with the other functionals can be found 
in the supplementary data. Small core pseudo-potentials 
were used for the heavy atoms (Mn and Y) associated 
with semi-valence and valence 2£ and 2£ plus polariza- 
tion basis sets 21 . The oxygen ions are represented in an 
all-electrons basis set of 2£ quality specifically optimized 
for O 2- ions^i. The mono-determinantal character of the 
Kohn-Sham approximation does not allow the possibility 
to deal with magnetic frustration as observed in YMnOs. 
Indeed, only collinear magnetism is possible in the CRYS- 
TAL code. A non-magnetic electronic density can be cal- 
culated, however it does not account for the open-shell 
and high spin character of the Mn 3+ ions, which are re- 
sponsible for a large stabilization energy (of the order of 
two Hubbard 3<i-shell U per Mn atom, this is of the order 
of lOeV) and crucial for a reasonable electronic structure 
representation. A "pseudo" antiferromagnetic configura- 
tion can be proposed on a 2 x 2 x 1 cell. In this case, only 
two among the three antiferromagnetic (AFM) couplings 
are satisfied within each triangle. It results in a lower en- 
ergy for the antiferromagnetic bonds and an associated 
loss of symmetry. When the geometry optimization is 
performed, the Mn-Mn distances associated with anti- 
aligned spins are shorter than the bonds with aligned 
spins, resulting in a structure with unphysical symme- 
try and characteristics. The choice of a ferromagnetic 



(FM) configuration for the Mn 3+ S=2 spins thus seems 
the best compromise. Indeed, not only this configuration 
respects the system symmetry (a crucial aspect for accu- 
rate phonons modes), but also the energy error that can 
be associated with the loss of the inter-atomic exchange 
energy (that is 1/2 J ~ 1.5 meV per bond) is much weaker 
than the error associated with the loss of the magnetic 
character for the Mn 3+ ions. For these reasons all our cal- 
culations will be done within a ferromagnetic Mn spins 
alignment. 

Within this procedure, we computed the polarization 
(using Berry phase approach) for the different experi- 
mental geometries available in the literature at different 
temperatures. Table [I] reports the computed polarization 
values. One sees immediately that our computed values 
are in very good agreement with the experimental ones, 
despite the fact that the magnetic order used in the cal- 
culation is ferromagnetic and not antiferromagnetic. The 
quality of these results thus validate our choices for the 
calculations. In addition it supports the idea, first pro- 
posed by Lee et afii, that in this system the magneto- 
electric coupling is mainly driven by elastic effects rather 
than by a direct coupling between the polarization and 
the magnetic order. 
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Table I: Spontaneous polarization Ps in \iC.cm~ 2 for different 
temperatures. Experimental and computed values using the 
experimental geometries at the given temperatures. 



In a second step we optimized the geometry. Both 
the cell parameters and the atomic positions compare 
very well with the experimental values in the paramag- 
netic phased 2 -. Indeed, the error on the lattice param- 
eters is weaker than 0.5% and for the atomic positions 

Eat K7 p -r a T lc | 2 = 0.0087 A 2 . 

Finally we computed the Born effective charge tensor 
Z* for the optimized geometry. They are reported in 
table [III Again our results compare well with the experi- 
mental evaluations validating our methodological choices. 



III. THE PHONONS SPECTRA 

Table Hill displays the T point optical transverse (TO) 
phonon modes, both for the paraelectric (columns 1-2) 
and ferroelectric (columns 6-7) phases. Raman and Infra- 
Red (IR) experimental data are given for comparison. 
The P^z/mmc group of the paraelectric phase decom- 
poses in 

lA lg © 3A 2u © 3B lg © 2B 2u © lE lg © 3E lu © 3E 2g © 2E 2u 

optical modes, out of which the A\ g , E\ g and E 2g are 
Raman active while the A 2u and E\ u are IR active, leav- 
ing the Bi g , B 2u and E 2u modes inactive both in Raman 
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Table II: Born effective charge. The dynamic charge, 
Tr(Z*)/3, is compared to the experimental evaluation from 
infra-red measurements 16 . 



and IR. The Pterin group of the ferroelectric phase de- 
composes in 

9Ai © hA 2 © 5£i © 10£ 2 © 14£i © 15£ 2 

optical modes, out of which the A±, E\ are active both 
in Raman and IR, the E 2 are active in Raman only and 
the A 2 , B\ and B 2 are inactive both in Raman and IR. 
For the ferroelectric phase, the optical longitudinal 
(LO) phonon modes have also been evaluated and are 
reported in the supplementary data. 



A. High temperature discussion 

At high temperature only four modes were experi- 
mentally measured by Fukumura et al 14 . These modes 
were respectively assigned by Fukumura et al to one A\ g 
mode (666 cm -1 ), one E\ g mode (420 cm -1 ) and two E 2g 
modes (120 and 395 cm -1 ) of the P6s/mmc group of the 
high temperature PE phase. In addition figure 5 of ref- 
erence LLJ clearly shows that these modes are observed 
continuously from 300 K to 1200 K in the same energy 
ranges. Going back to our calculations, if we want to 
assign these modes by continuity from the low tempera- 
ture phase (P6scm to P6s/mmc)^ we have to suppose 
that the mode at 395 cm -1 should be assigned as E 2u 
(computed at 402 cm -1 in the P6s/mmc group and at 
392 cm -1 in the P6scm group, and measured at 406 cm -1 
in the low temperature phase) and not as E 2g (computed 
at 500 cm -1 in the P6s/mmc group and at 473 cm -1 in 



the P6scm group, and measured at 483 cm -1 in the low 
temperature phase). Similarly, the mode at 420 cm -1 
should be assigned as E\ u and not as E\ g and the mode 
at 120 cm -1 should be assigned as E 2u rather than E 2g . 
At this point there is a clear contradiction between the 
continuity requirement of the phonons modes between 
the low and high temperature phase and the Raman sym- 
metry requirements. Indeed the "ungerade" E\ u and E 2u 
irreducible representations (irreps) are not Raman active 
and no crystal disorientation or twinning could induce an 
inversion between a "gerade" and an "ungerade" mode. 
We must thus assume that the symmetry group at 1200 K 
cannot be the P6s/mmc group. 

What type of hypotheses are left? The first possibility 
is the existence of an intermediate phase in the 1200 K 
temperature range as suggested by several authors ^ 25 i 26 . 
However the PQ^/mcm symmetry group suggested is no 
more than the P^/mmc group, compatible with the 
phonons spectrum. Indeed, symmetry group analysis 
shows that the P6s/mmc to P6s/mcm transformation 
only exchange the B\ g with the B 2gi and the B\ u with 
the B 2u irreps. A second possibility would be that the 
intermediate phase corresponds to the P63171C group. In 
this group the inversion center is lost, and the E\ g and 
E\ u irreps of P6s/mmc are associated with the E\ ir- 
rep of the P63171C group and with the E\ irrep of the 
low temperature phase (resp. E 2g and E 2u are associ- 
ated with E 2 ). Both the E\ and E 2 irreps are Raman 
active in the P63171C group and the experimental Fuku- 
mura et a'— data can be easily associated with the com- 
puted modes (see table IIIIj, columns 3-4). The average 
error between the computed and measured frequencies 
is 18 cm -1 , essentially supported by the higher A mode 
(3 cm -1 for all the other modes). The main problem with 
this hypothesis is that the unit cell tripling at the P6smc 
to P6scm transition should have been seen in diffraction 
experiments, while it does not seem to be the case^. 
Therefore, the only possibility compatible with all data 
seems to be that a P6s/mmc to P6scm phase transition 
occurs at a temperature higher than 1200 K, as proposed 
by Jeong2£ and Gibbs^ from neutrons diffraction mea- 
surements (T c = 1258 K ± 14 K). In other terms, all the 
high temperature data from Fukumura et al 14 should be 
interpreted within the ferroelectric P63C171 group. 



B. Low temperature discussion 

At low temperature our computed modes fit quite well 
with the available experimental data, both Raman and 
IR (see table HTH columns 6 and further). Table [TV] dis- 
plays the average errors between our computed modes 
and the experimental data for each irreducible represen- 
tation. The experimental phonons modes can be associ- 
ated with the computed ones except for a few exceptions. 
i) The mode at 190 cm -1 seen by Iliev et al in the A\ irre- 
ducible representation cannot be associated with a com- 
puted mode. This mode, not seen by the other authors, 
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Table III: Hexagonal YMnOs vibrational frequencies in cm -1 for the paraelectric (P6s/mmc) and ferroelectric (P6scra) phases. 
The frequencies for one of the possible intermediate phases (PQzmc) was also computed. Modes with problematic assignment 
are in bold faces ; modes that are symmetry allowed but not seen are designed by dashes. 
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Table IV: Average errors per irrep on the YMnOs vibra- 
tional frequencies compared to different experimental data 
(in cm -1 ). 



cannot be associated with a A\ mode. Indeed, the A\ 
modes are IR active and all predicted modes are observed 
by Zaghrioui et ai^ up to more than 450 wave numbers. 
If not artifactual, the 190 cm -1 mode must belong to an- 
other irreducible representation. The only possibility in 
this energy range is an E2 mode predicted at 150 cm -1 . 
Indeed, all other IR or Raman active modes in this energy 
range are experimentally observed and easily associated 
with the computed phonons. ii) The 210 cm -1 A\ mode 
seen by Fukumura et aS^. Once again a A\ mode around 
210 cm -1 is incompatible both with the computed values 
and with the other experimental results (IR or Raman). 
Most probably this mode has been improperly assigned 
as a A\ mode and should rather be considered as a E\ 
(seen at 211 cm -1 by Zaghrioui et aft£) or a E2 mode 
(seen at 215 cm -1 by Iliev et a&£). 

Another point we would like to address is the reason 
why all phonons modes are not observed in experimental 
data. As far as the IR active modes are concerned we 
computed an estimate of the intensities, the analysis of 
which gives us some insight on the reasons some modes 
are not observed. In the A\ representation two modes are 
not seen in reference [l6|, namely the highest mode com- 
puted around 720 cm -1 and the mode computed around 
460 cm -1 . The IR computed intensity of the 720 cm -1 A\ 
mode is very weak (between 5 and 12 km/mol according 
to the functional used) and it is thus not surprising this 
mode is not seen experimentally. The IR intensity of 
the A\ mode computed around 460 cm -1 mode is of the 
same order of magnitude as the next computed A\ mode 
at 492 cm -1 which was associated to the experimental 
mode found at 434 cm -1 ). Even if the IR intensity re- 
mains small it is non negligible and the second mode is 
observed. Assuming an equivalent shift in energy be- 
tween experiments and calculation on these two modes, 
the first mode should be searched around 400 cm -1 in 
the IR spectrum however it is not seen. In the E\ ir- 
reducible representation, five modes are not seen in IR 
experiments, namely the modes computed at 269 cm -1 , 
386 cm -1 , 461 cm -1 , 515 cm -1 and 668 cm -1 . The calcu- 
lations reveals that all these modes are associated with 
very low IR intensities that could explain there are not 
seen by Zaghrioui et a££. 



C. The A2 phonon mode at 33 cm x 

We would now like to discuss in a little more details 
the first optical mode, computed at 33 cm -1 ~ 4meV 
in the A2 irreducible representation. The displacements 
vector associated with this phonon mode is reported on 
table [V] and pictured in figure [2j 
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Table V: Displacements vector (normalized to classical am- 
plitudes) associated with the first A2 mode found around 
33 cm -1 . 




Figure 2: (Color online) Schematic picture of the atomic 
displacements in the (a, b) plane for the lowest Ai phonon 
mode. The blue, solid lines, Mn triangles remain nearly un- 
touched and rotate as a whole around the c axis. The red, 
dashed and magenta, dotted Mn triangles are alternatively 
enlarged/reduced, the Mn ions moving within the (a, b) plane 
while the central oxygens move conversely along the c di- 
rection. A jmol xyz file can be found in the supplementary 
data in order to visualize the displacements of one manganese 
plane. 

This phonon mode, active neither in Infra- Red nor in 
Raman experiments exhibits displacements that can be 
expected to strongly affect the magnetic exchange cou- 
pling constants along the red, dashed and magenta, dot- 
ted bonds as picture in figure [2j 

Let us notice that the effective exchange J between 
two manganese atoms is the sum of a direct exchange 
contribution Jj (Pauli exchange, always ferromagnetic, 



exponentially dependent on the Mn-Mn distance) and 
a through-ligand super-exchange term J\ (antiferromag- 
netic, dependent on the metal ligand distances and on the 
Mn-O-Mn angles). In YMn03 the magnetic exchange 
couplings are known to be antiferromagnetic showing 
that the through-ligand super-exchange term J\ domi- 
nates over the direct exchange contribution J^. We ana- 
lyzed this super-exchange term and reported in figure [3] 
the main super-exchange paths. J\ thus scales as the 
sum of the four associated super-exchange terms, each of 
them scaling as the square of the product of the associ- 
ated Mni3d — 02p orbital overlap and the 02p — Mn 2 3d 
orbital overlap. 






(Jh) 



(Ji 



(j h 






Figure 3: Main super-exchange paths between two manganese 
atoms. Let us remember that in YMnOs the manganese is 3d 4 
and the singly occupied orbitals are d xy , d x 2_ y 2, d xz , d yz if z 
is along the c axis. We can evaluate 

J h ~ {Mmd xz \Op z ) 2 {Op z \Mn 2 d xz ) 2 
J h ~ (Mmd yz \0 Pz ) 2 (0 Pz \Mn 2 d yz ) 2 
J h ~ {Mmd x 2_ y 2\Opy) 2 (0 P y\Mn 2 d x 2_ y 2) 2 
Ji A ~ {Mmd xy \Op x ) 2 (Op x \Mn 2 d xy ) 2 

The A2 phonon mode under consideration involves the 
following atomic displacements on the red, dashed and 
magenta, dotted Mn triangles (see figure [2]). 

• The Mn-Mn bond length elongation / contraction 
(in opposite phases between red and magenta trian- 
gles). This movement affects the direct part of the 
exchange : J^, increasing its ferromagnetic charac- 
ter when the Mn-Mn distances are contracted and 
decreasing it when the Mn-Mn distances are elon- 
gated. 

• The up and down movement (along the c direc- 
tion) of the oxygens atoms at the center of these 
triangles. This movement strongly affects both the 
Mn-0 distances, reducing them when the oxygens 
come closer to the Mn plane, and the Mn-O-Mn 
angle. 

One sees immediately using figure [3] that, when the oxy- 
gen comes closer to the Mn plane, the overlap between 



the magnetic orbitals of the Mn atoms and the 2p bridg- 
ing orbitals of the central oxygen increases, thus increas- 
ing the antiferromagnetic character of J\. At the same 
time the Mn-Mn bond length increases thus decreas- 
ing the ferromagnetic character of Jd- Both terms goes 
in the same direction resulting in a enlarged antiferro- 
magnetic character for the bonds of the (alternatively 
red/magenta) triangles when the central oxygen comes 
closer to the Mn plane and a reduced antiferromagnetic 
character when the central oxygen goes away from the 
Mn plane. 

The blue, solid line, Mn triangles (centered around the 
unit cell origin) are subject to a rotation as a whole 
around the c axis going through their central oxygen 
(O3). One thus does not expect the magnetic exchange 
within the blue Mn triangles to be affected by this A2 
phonon mode. However this rotation results in a strong 
distortion of the triangles with blue red and magenta 
bond inducing a magnetic interaction scheme quite far 
from the nearly homogeneous triangular lattice of the 
static structure. 

The main conclusion of this analysis is that this A^ 
phonon mode is strongly coupled to the magnetic inter- 
actions. 

Let us now put these results into perspective with the 
characteristics of the hybrid phonon-magnon excitation, 
studied using inelastic polarized and non-polarized neu- 
tron scattering by Pailhes and coworkers 4 . One notices 
that this A2 mode seems to fulfill all the requirements 
for a good candidate for the phonon participating to the 
hybrid mode. Indeed, it is an optical mode in the ade- 
quate energy range, exhibiting atomic displacements con- 
sistent with the umbrella motion proposed in reference |4J 
and that can be expected from the above analysis to be 
strongly coupled to the spin degree of freedom. 



IV. CONCLUSION 

The present paper proposes T point phonon calcula- 
tions for the YMnOs compound both in the ferroelectric 
phase and in the paraelectric phase. Our calculations 
agrees well with the Infra-Red and Raman experimental 
data. The experimental versus calculated modes corre- 
spondence is discussed in details for the few problem- 
atic modes. We were able to explain the fact that sev- 
eral modes could not be seen in IR experiments due to 
their very low (calculated) intensities. A careful anal- 
ysis of the phonons modes correspondence, when going 
from the ferroelectric to the paraelectric phase, is per- 
formed and leads to the conclusion that, the phonon 
modes observed in Raman scattering at 1200K 14 can- 
not be associated with the paraelectric phase P6s/mmc 
space group. The different intermediate subgroups al- 
lowed between the paraelectric P6s/mmc space group 
and the ferroelectric P6scm group were also checked 
against the modes continuity and the experimental data 
and had to be discarded. The only possibility agreeing 



with the phonon calculations and the experimental data 
is that the P^/mmc to PQ^cm phase transition occurs at 
a temperature higher than 1200 K, supporting the propo- 
sition of Jeong 27 and Gibbs 8 (T c = 1258 K ± 14 K). 

Finally we studied in more details the first optical 
mode. This mode belongs to A2 irreducible representa- 
tion and is inactive both in IR and Raman experiment. 
This mode however seems to have all necessary charac- 
teristics to be a good candidate for the phonon part of 
the electro-magnon mode observed in inelastic neutrons 
scattering by Petit, Pailhes et a&£. 
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